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Abstract 



We show that an instabihty may be present in resonant tunnehng through a 
quantum well in one, two and three dimensions, when the resonance lies near 
the emitter Fermi level. A simple semiclassical model which simulates the 
resonance and the projected density of states by a nonlinear conductor, the 
Coulomb barrier by a capacitance, and the time evolution by an iterated map, 
is used. The model reproduces the observed hysteresis in such devices, and 
exhibits a series of bifurcations leading to fast chaotic current fluctuations. 
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Resonant tunneling through a double barrier, originally a rather elementary academic 



and surprising physics.!!! One instance is the bistable region in the I-V characteristic that 
some samples exhibit, leading to hysteresis as the curve is traced first increasing and then 
decreasing the bias.oi This effect is currently understood as caused by the interaction of 
the current flowing through the device with the charge trapped in the well formed between 
the barriers.! Calling Q this charge, its effect on the incoming electrons may be viewed as an 
increase in the local potential by Q/C, where C is the capacitance, lifting the resonance level 
in the same amount. The latter can thus be somewhat above or below the emitter conduction 
band edge at one same bias, depending on whether the well is charged or uncharged. Since 
in the first case there is current flow through the well and in the second there is none, both 
cases are physically consistent, and in fact, observed experimentally. As an application, it 
has been sioeeested that in the bistable region the device may act as a THz detector and 
fast switch.B'ty 

Another instance are the THz oscillations that may arise in the presence of a magnetic 
field in the direction of the current flow, which become chaotic if the field is sufficiently 
strong.lilhEi In contrast with the hysteresis effect described above, these oscillations are 
associated with passage of the resonance across the emitter Fermi level. Assuming the 
resonance is initially above the emitter Fermi sea and off resonance, it will eventually enter 
the latter as the bias is increased, allowing electrons to tunnel into the well. Thus, with 
a certain time constant, charge begins to build up between the barriers. The Coulomb 
repulsion between this charge layer and the incoming electrons effectively cause an upward 
shift of the bottom of the well, and hence, also of the resonance state. This may drive 
the system off-resonance again until the charge, due to the finiteness of the barrier and the 
associated lifetime of the resonance state, has tunneled out again into unoccupied states 
on both the emitter and collector sides. The resonance then is no longer sustained above 
the Fermi energy by the decreased charge in the well, the resonance falls again and a new 
cycle begins. Whether this oscillation is damped away or not is determined by the strength 
of the interaction, which in turn depends on the projected density of states in the emitter 
at or near the Fermi energy. By changing the profile in the density of states the magnetic 
field enhances this coupling, and, if large enough, induces the system to perform sustained 
oscillations. It is the purpose of this work to show that, under favorable conditions, the 
oscillations may also be present in the absence of a magnetic field, as well as situations of 
lower dimensionality such as tunneling through a quantum wire or dot. 

We model the system by a simple circuit that captures the essence of the system under 
discussion, shown in Fig. 1. As will follow from the discussion below, it includes the following 
features that appear to be essential in the description of the device: (1) the existence of a 
resonance state with an associated lifetime, which results in a time delay between the bias 
voltage and the build-up or decay of the charge in the well ("tunneling time"); (2) the 
existence of a Fermi sea on the emitter side, with a density of states depending on the 
dimensionality of the device (bulk, planar or linear), and on the magnetic field, if present; 
(3) the existence of the Coulomb repulsion between the charge trapped in the well and the 
incoming electrons. The nonlinear element in the circuit represents transmission from the 
emitter side into the well, allowing a charge 
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to leave the source Uq in the time interval r. The dimensionless nonlinear function f{V), 
depending on the voltage V = Uq — U between emitter and well, is the convolution of the 
transmittivity of the resonance with the number of occupied states in the emitter, available 
for tunneling. This latter quantity may be written as 

N{V) = e{eV -Er + Ef) [1 - e{eV - Er)] Pd{V), (2) 

where 9{x) is the Heaviside step function, Ej. is the energy of the resonance level at zero 
bias and Ep, the Fermi energy. All energies and voltages are assumed to be in electron-volts 
and measured with respect to the bottom of the conduction band on the emitter side. The 
quantity PdiY) represents the number of states in the emitter side, with the component of 
kinetic energy along the current flow equal the resonance energy shifted by the bias V. It is 
given by 

Pd{V)^aa{eV + EF-Er)'-^. (3) 

Here aa is a constant depending on the dimensionality ci = 3, 2 or 1 of the emitter (the well 
has dimension d — 1). The values are, cti — 2,a2 — CKs = m* L"^ / nh'^ , where L 

is the width of the emitter and m* the effective mass of the carriers. 

The capacitance C represents the Coulomb barrier. As charge flows through the nonlinear 
element, entering the well, the voltage drop V = Uq — U is reduced. At the same time charge 
is allowed to leave the well through the collector represented by the load resistor R. These 
elements define a time constant tq — RC characterizing the rate at which the well may be 
emptied. With these definitions the conservation of charge provides an equation of motion 
for the voltage U (t) on the output side of the device, 

C{U{t + T)-U{t)) = AQ-T^. (4) 
Defining Vn = V{t + nr) one obtains from Eqs. (1) and (4) the following iterated map, 

K+l = (1 - 7)K - iRImaJiVn) + lUo, (5) 

where j = t/tq and we assume that Tq > t holds. The map has fixed points given by 

V*^Uo-RImaxf{V*). (6) 

A simple linear analysis in the neighborhood of one of these fixed points shows that for it 
to be stable the condition 

-l<i?7^„,^|y, < 2/7-1 (7) 

must be satisfied. It follows that for either an increasing or decreasing function f{V) at 
the fixed point, an instability occurs provided the local derivative is sufficiently large in 
absolute value for the above condition to be violated. For a sharp resonance the function 
f{V) will essentially follow the profile given by Eq.(3) within the window E^ — Ep <V<Er, 
becoming negligibly small elsewhere. The finite resonance width smoothens the edges defined 
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by the top and bottom of the Fermi sea. Assuming a lorentzian form for the transmissivity 
T{V) = rV[(E -Er + eVf + r^] one has for V « Ep^ 

fiV) ^ p,iV) [arctgi ^) - arctg{ ^ . (8) 

The expression in parenthesis is a positive definite funtion of V resembhng a hat, giving f{V) 
such shape for d=l since then Pd{V) is a constant; for d=2 and 3 f{V) rather resembles an 
asymmetric hat. For convenience we call it the "hat function" in what follows. 

The fixed points, as determined by Eq. (6), may be found graphically as the intersection 
of the straight line through the origin V*, and the curve representing the right hand side. 
Because of the inverted hat form of the latter one can easily verify that there are either one 
(for large or small value of f/o), or three, fixed points (some intermediate region). In this 
last case the fixed point in the middle intersects the rather steep fall of the arctg function, 
yielding a likely unstable fixed point. The intersection furthest to the right (RFP) occurs 
in the flat portion of the hat function and is therefore stable. 

We are interested in the leftmost fixed point (LFP) of the triple solution case, representing 
the intersection located in the left side of the hat function, which may also be present when 
there is just one fixed point. As f/o varies, the intersection closely follows the functional 
form of the number of states available for tunneling (3). We first consider the d = 1 case, 
for which pi{V) is a constant in the relevant region. Then the LFP traces the contour of 
the arctg function and again the fall is governed by the rather steep arctg, and, as with the 
RFP, will generally violate the stability condition (7). The situation is as in the d=3 case 
with magnetic field since the latter modifies the spectra in the emitter transforming it into 
a series of quasi one dimensional dispersion laws, and for which the instability is known to 
exist .0 For d = 2, there is a square root dependence at the intersect, and if F is sufficiently 
small the diverging local derivative will cause the stability condition again to be violated 
near the edge (the Fermi energy). 

The d = 3 case needs closer attention because of the linear form of psiV). Calculating 
the total emitter current flowing through the resonance and relating the time constant r 
with the resonance width through r ~ 2TTh/T one gets, away from the hat edges 

df e^7Ti*D 

Rlmax^W* - (9) 

where Dz is the separation between the barriers, and e, eo are the dielectric constants of 
the enclosed material and the vacuum, respectively. Using the values appropriate for GaAs 
(m* = 0.068m, e = 12.5) one gets the simple stability condition 

1.3D^(nm) < 2 -7, (10) 

where the distance is to be given in nanometers. This condition is easily violated in 
actual samples. 

We have done numerical studies in order to verify the presence of the instability, and 
have found it to occur in all three cases discussed above. Figure 2 shows the I-V curves for 
the (a) d=l, (b) d=2 and (c) d=3 cases in the absence of a magnetic field. The parameters 
used in the figure are 7 = 0.5, F = ImeV, Ep = 20meV, E^ = ASmeV and = 13nm, 
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appropriate for typical devices based on AlGaAs. It is also necessary to specify some short 
cross sectional dimensions in cases (a) and (b), chosen as D^Dy = lOOnm^ and = lOnm 
in the figure, respectively. For these values, at low bias Uq the device is always stable. Yet, 
as the voltage increases and current starts to fiow, marking the entrance of the resonance in 
the Fermi sea, the iteration does not settle to a fixed point but rather oscillates in steps r, 
first regularly (bifurcations region) and then, at larger bias, in a chaotic fashion. The case 
with magnetic field and d = 3, not shown, produces a similar figure for large enough field, 
as expected from results reported previously.0lll 

In summary, we have shown that an instability induced by the interaction of the current 
with the trapped charge may appear in resonant tunneling through constrictions in one, 
two and three dimensions, when the resonance is close to the Fermi level. Instead of inte- 
grating the quantum mechanical equations of motion as done previously for the d=3 case 
with magnetic fieldS, we have used a circuit model with the advantage of mathematical 
simplicity and the possibility of an analytical discussion of the stability in one, two and 
three dimensions. Experimental accesibility of the instability may require testing emission 
or absorption of THz radiation. Usual electronics will normally just register a time average 
of the oscillations because of the small value of their typical period, which, for a resonance 
of ImeV width in a GaAs quantum well would be of the order of 4ps. 

This research was supported in part by a Catedra Presidencial en Ciencia and FONDE- 
CYT grant 1990425. 
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FIGURES 



FIG. 1. Model circuit for a double barrier device. The function f{V) represents the response 
for tunneling into the well 

FIG. 2. Current-voltage iterates for a quantum well embedded in a (a) one dimensional, (b) 
two dimensional and (c) three dimensional conductor 
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